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We present a microscopic approach to dynamical pionisa- 
tion of the hot quark-gluon matter formed in ultrarelativistic 
heavy ion collisions. The time evolution of the system is de- 
scribed assuming that quarks undergo Brownian motion in a 
thermal bath provided by the gluons. The rate of hadroniza- 
tion as well as the time dependence of the temperature of the 
system are seen to be quite sensitive to the QCD A param- 
eter. Even in a non-equilibrium scenario, we find that there 
appears a clear hint of a first order phase transition. 
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Search for the novel state of strongly interacting matter 
called Quark-Gluon Plasma (QGP), in which the prop- 
erties of bulk matter are determined by the quarks and 
gluons rather than the usual hadrons, is a major sphere of 
current activity. Quantum Chromodynamics (QCD) pre- 
dicts that such a phase, i.e., a metastable configuration, 
should exist at very high temperature and/or baryon den- 
sity. It is strongly believed that conditions conducive to 
the formation of QGP could indeed be achieved, albeit 
transiently, in highly energetic collisions of heavy ions 
and consequently, a substantial effort is now being de- 
voted to look for QGP in such collisions. 

An inescapable feature of the collision process is that, 
because of colour confinement, the quarks and gluons 
in the QGP must, at some epoch, turn into hadrons 
which would ultimately be detected, never the individ- 
ual quarks and gluons. An ab initio calculation of this 
very important, non-perturbative, phenomenon is, un- 
fortunately, beyond the capability of QCD as yet. It 
has thus been widely postulated that this conversion of 
coloured particles to the colourless hadrons en bulk could 
be an actual phase transition (the order of which is an 
open issue) |Q. Assuming the phase transition to be of 
first order, some authors Q have attempted to investi- 
gate the process of hadronization in a QGP through a 
bubble nucleation mechanism. In addition to being a 
classical picture, this mechanism also requires thermo- 
dynamic equilibrium. Recent studies on QGP dynamics 
j^-H] , showing the lack of thermodynamic equilibrium in 
the quark-gluon phase in ultrarelativistic heavy ion col- 
lisions, indicate that such an ideal situation is rather un- 
likely. 



Recently there have been some attempts to study the 
formation of hadrons in quark matter using different 
semi-microscopic approaches. These studies can be char- 
acterized either as model dependent calculations [||, or 
the computer codes based on the string phenomenology 
, or other phenomenological description of hadroniza- 
tion J^] . None of these approaches account for the essen- 
tial lack of equilibrium in the quark-gluon phase. Some 
efforts have also been made to estimate hadronization 
within the parton cascade model by introducing a cut-off 
to mimic the non-perturbative effects To the best of 
our knowledge, the first study aimed at investigating the 
dynamical process of hadronization in a non-equilibrated 
quark-gluon system from a physically transparent ap- 
proach was in In this letter, we present the de- 
tails of this scenario describing the evolution of the non- 
equilibrated hot quark-gluon matter towards hadroniza- 
tion. For simplicity, we have restricted ourselves to the 
case of pions only in the present work; since the pions 
account for the bulk of the multiplicity, the conclusions 
should remain valid not only qualitatively and even semi- 
quantitatively. 

We have earlier shown that perturbative estimates 
of the gluon-gluon, quark-gluon and quark-quark cross 
sections allow us to study the evolution of the quark- 
gluon matter formed in ultrarelativistic heavy ion col- 
lisions by visualizing the quarks as Brownian particles 
in a hot gluonic thermal bath. In this work, we start 
with such a premise, which, in our opinion, describes the 
non-equilibrium aspects of the evolution in a physically 
transparent manner. 

Hadronization in such a system can be studied in the 
light of Smoluchowski's theory of coagulation in colloids, 
which was elaborated further by Chandrasekhar [ jTlJ . 
This theory suggests that coagulation results as a con- 
sequence of each colloidal particle being surrounded by 
a sphere of influence of a certain radius R such that the 
Brownian motion of a particle proceeds unaffected only 
so long as no other particle comes within its sphere of 
influence. When another Brownian particle does come 
within a distance R of the test particle, they form a two 
body cluster. This cluster also describes a Brownian mo- 
tion but at a reduced rate due to its increased size/mass. 
The process continues further till a single cluster of all 
the particles is formed. 
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In order for this stochastic scenario of cluster forma- 
tion to apply to a system of Brownian quarks in the hot 
gluon bath, it is essential that each quark has an appro- 
priate sphere of influence of radius r. Obviously, this 
radius r will depend on the spin-isospin combination of 
final cluster (whether the final cluster is scalar, pseu- 
doscalar, vector or axial vector meson or even a baryon). 
Mesons are formed when one quark and one antiquark 
with proper quantum numbers come within the spheres 
of influence of each other. (Clusters with greater num- 
bers of quarks and antiquarks (e.g. baryons) can also 
be formed by imposing the conditions of colour neutral- 
ity and charge balance properly.) This implies that the 
radius of the sphere of influence corresponds to the corre- 
lation length between the quarks in the proper hadronic 
channel. In other words, this is the screening length of 
the corresponding hadrons in the hot quark-gluon mat- 
ter. For the present purpose, the radius corresponds to 
the screening length of pions, which can be evaluated in, 
say, the Nambu-Jona-Lasinio (NJL) model |12| or the lat- 
tice |l3|| . (One can immediately see, without further ado, 
that the rate of pion formation in the hot quark matter 
should be rather small at high temperatures and increase 
with falling temperature as the pion screening length (in- 
verse of the screening mass) decreases with increasing 
temperature in all realistic pictures.) The pions formed 
at very high temperatures are most likely to decay back 
into quarks and antiquarks. It has been recently shown, 
within the NJL model, that the pion decay width is very 
large at high temperatures and becomes zero around a 
temperature of 215 MeV or so |L4|]. Since this would 
be a crucial ingredient in what follows, we adopt, for the 
sake of consistency, the NJL model estimates for the pion 
screening mass in the present work. Note that evaluating 
the pion decay width as a function of temperature on the 
lattice is not possible as the dynamical pion mass cannot 
be evaluated there. 

The rate of pion formation from Brownian quarks as a 
stochastic process, as is evident from the preceding dis- 
cussion, depends on the number of quarks (antiquarks) 
falling into the sphere of influence of another antiquark 
(quark), the number of pions decaying back to quarks and 
antiquarks and also the change in the pion density due 
to the expansion of the system. Thus the rate equations 
can be written as, 
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In eqs.(|i],|2|), the first term is the rate of pion formation 
(a = + or -); the second term is the rate of pions decay- 
ing back to quarks and the third term is due to Bjorken 
(longitudinal) expansion of the system. i(J) stands for u 
or d (we ignore s and other heavier flavours). < v ■ f > 
(the average relative velocity in the radial direction) is 
calculated using the Jiittner distribution, 
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There would also be a corresponding rate equation for 
the antiquarks, which looks exactly like eq. (0) and hence 



not explicitly written. In eq. the T 



total 



as well as 



Y total gg^qq stand for the corresponding net quantities. 

As mentioned earlier, we are considering a non- 
equilibrated quark matter and hence the pions formed 
will also be out of equilibrium. This is taken into ac- 
count by multiplying the relevant distribution functions 
with the ratios r q = n q /n eq and rv = n^/n^ where n q 
and n eq are non-equilibrium and equilibrium densities of 
quarks and and n e7r are non-equilibrium and equilib- 
rium densities for pions. These details are given in 

In all these expressions, the appropriate masses are the 
effective masses including the current as well as thermal 
contribution, , whose importance in determining the dy- 
namics of the hot quark matter has been well established. 
For quarks (antiquarks), this is 
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and m q (curr) is taken to be 10 MeV. For gluons the 
thermal mass is, 



m g (thermal) — —g s T 
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The running coupling constant a s as a function of tem- 
perature is given by |0 



12tt 



(33 - 2n f )ln 



(7) 



with Q 2 = m eff 2 (T) + 9T 2 . 

Simultaneously, we must take account of energy mo- 
mentum conservation which, for a Bjorken flow, corre- 
sponds to the following equation 
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where e = e to 



e n . We also include the one 



loop correction to e g MM . e and P are related through the 



2 



velocity of sound, as in H. For a complete description 
of the system, eqs. (P, (0), (||) and (||) must be solved 
self-consistently. The initial conditions are taken from 
|^| for RHIC energies. The initial time (t g ) is the time 
when gluons thermalise (=0.3 fm), where r q — 0.15, r g =l 
and r n is taken to be 0. The temperature at this time 
is 500 MeV. The pion decay width, dynamical mass and 
screening mass are taken, as already mentioned, from the 
NJL model ]l2||l4| ]. Note that we are working at y = 
so that t and r are the same and the baryon chemical 
potential is zero. 



Figure 2 shows the variation of temperature with time. 
Obviously, there is a dramatic effect of the QCD param- 
eter A. In all the cases, there is a change at T ~ 215 
MeV, corresponding to t ~ 3.5 fm; the variation of tem- 
perature with time becomes slower, as is expected in the 
mixed phase of a first order phase transition. At A=0.2 
GeV, this occurs for a very short period of time, before 
the system starts cooling again. The duration of the con- 
stant temperature configuration increases with A, and for 
A=0.4 GeV, it persist upto 9 fm before the temperature 
of the system starts falling again. 
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Figure 1 : Time evolution of quark and pion number 
densities for various values of QCD parameter A. 

Figurel shows the variation of pion and quark num- 
ber densities with time, for various values of the QCD 
parameter A. In all three cases, we find the same qual- 
itative feature that pions start appearing in the system 
quite early on but they become appreciable in number 
only after some time. At late times the system is domi- 
nated by pions. This cross over occurs at t > 4 fm for A 
= 0.2 or 0.3 GeV while for A= 0.4 GeV this happens at 
t ~ 6 fm. 
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Figure 2 : Time evolution of the Temperature for 
various values of QCD parameter A. 

Obviously, this is a clear indication of an apparent first 
order transition. Microscopically, the appearance of the 
mixed phase at a temperature of ~ 215 MeV can be un- 
derstood from the fact that the pion decay width goes 
to zero at such a temperature All the pions that 

were formed earlier in the system tended to decay back 
to quarks and antiquarks on a fast time scale. Only after 
the pion decay width becomes small would the formed 
pions become stable. 

The role of the QCD parameter A is better understood 
from figure 2. The higher the value of A, the higher the 
quark thermal masses. As a result, the lower the relative 
velocity, which would lead to a lower rate of pion produc- 
tion as well as a slower depletion of the quark number 
density. Thus the mixed phase, i.e. the domain where 
quarks and pions have comparable densities, would not 
only occur later in time but also persist for a longer du- 
ration. It should however be mentioned that there are 
also other effects associated with the QCD A (like the 
quark production rate from gluon fusion and /or decay, 
the gluon energy density and so on) which have com- 
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peting roles in determining the actual number densitites. 
This may be the reason why one does not notice a drastic 
difference between A = 0.2 and 0.3 GeV in figure 1, while 
at A = 0.4 GeV, the change is more noticeable even in 
figure 1. A detailed analysis of all these different effects 
is in progress now. 

We must stress the fact that the first order phase tran- 
sition is only an approximate one. There is no unique 
demarcation between the different phases. High temper- 
ature domains are dominated by quarks with a few pions 
present in the system, the intermediate region (T ~ 215 
MeV) having comparable numbers of quarks and pions 
and the low temperature domain being mostly pions with 
a few quarks. Nonetheless, it can be fairly concluded 
from these results that the concept of a first order con- 
fining phase transition is not a bad approximation for 
the process of hadronization in QGP. It should also be 
noted that the so called critical temperature of ~ 215 
MeV, which corresponds to the temperature where the 
pion decay width vanishes, is a model dependent quan- 
tity. In our case this derives from the NJL model, but 
the qualitative features should of course be model inde- 
pendent. 

In conclusion, we have studied, in a physically trans- 
parent picture, the dynamical process of hadronization 
(pion formation) in a non-equilibrated quark-gluon sys- 
tem formed in ultrarelativistic heavy ion collision. Our 
results show that even in a microscopic analysis, there is 
a clear indication of an apparent first order phase transi- 
tion in the system. Nonetheless, the persistence of non- 
perturbative modes in the high temperature phase seems 
to be a real possibility, as also the existence of some 
colour degrees of freedom till lower temperatures. These 
issues deserve urgent attention in the context of QGP 
diagnostics. 
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